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ABSTRACT 

This paper 1 examines the relationship between the radar 
backscatterlng coefficient of a vegetation canopy, cr° can , and the 
canopy’s leaf-area Index (LAI). The relationship fs established through 
the development of a model for corn and sorghum and another for wheat. 
Both models are extensions of the cloud model of Attema and Ulaby 
(1978). Analysis of experimental data measured at 8.6, 13.0, 17.0, and 
35.6 GHz Indicates that most of the temporal variations of c^an can be 
accounted for through variations In green LAI alone. If the latter Is 
greater +han 0.5. 
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1.0 -LNIRQDilCIlQN 

Among the prime objectives of agricultural remote sensing Is the 
early and accurate estimation of agricultural production In a manner 
superior to that of more conventional means. The first step In that 
direction Is the Identification of crop types, which usually is 

accomplished by means of multi-date observations of the particular test 

* 

site or area of Interest. Once the Identities of the various 
field-covers In a given area have been determined, crop acreage can be 
estimated. In order to estimate crop production reliably, an estimate 
of ySe I d-per- unit-area Is needed for each of the crop types or fields 
concerned. In Section 2, traditional as well as remote- sens I ng 
approaches to crop yield estimation are discussed. 

Several papers have been pub I I shed over the past few years 
documenting the ability of radar to Identify crop types correctly, I.e., 
with a classification accuracy of 90 percent or higher, based on two or 
three multi-date observations (Bush and Ulaby, 1978; LI etal., 1980; 
Hoogeboom et a|., 1982; Brisco and Protz, 1981; and Shanmugan et al., 
1983). The purpose of this paper Is to examine the 
relationship oetween the radar backscatter Ing coefficient of a 
vegetation canopy ( a can ) and the canopy’s leaf-area- Index (LAI), the 
latter being related to the solar radiation Intercepted by the canopy, 
which In turn Is a fundamental component of crop-yield models. The 
means of such an examination Is based on an analysis of experimental 
data In terms 'fa modified form of the canopy "cloud" model of Attema 
and Ulaby (1978). 
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2.0 jSBQEJLLELCL . AND..LA I 

The yield of a particular crop Is a function of many variables. 

Including the nutrients available to It and the weather conditions 

affecting It over the growing season. Successfully estimating a crop’s 

yield requires knowledge of factors that Include the health of the 

plants and their vigor at various points Jn time, the amount of 

fertilizer and water available to them, and hourly or dally 

weather-condition Information (Daughtry etal., 1980). Me+eorol oglcal 

Information Is available from organizations such as NOAA for the United 

States and related regions, and from the World Meteorological 

* 

Organization (WMO) for worldwide weather monitoring. 

Early attempts to predict yields were based upon such factors as 
previous years’ yield, Improvements In technology, and the effects of 
current weather conditions (MacDonald and Hall, 1980). By Incorporatl ng 
A pr lor I knowledge such as soil type and water-table level, estimates of 
soil moisture Information are added to the model, thus improving Its 
predictive capacity. The amount of available solar radiation, which Is 
the basic source of energy In the process of photosynthesis, Is also a 
major determinant of overall yield. 

As an energy source, solar radiation Is available to tl.j plant only 
when It Interacts with the leaves. Thus, an estimate of the total area 
of exposed and active leaves, In conjunction with Incident solar 
radiation, provides another measure of crop performance, and 
consequently of final yield. This was the reasoning used by Dale (1977) 
In the development of the energy-crop-growth (ECG) model. This model 
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takes the form: ^ POOR QUAU'tY 

t Z 

ECG “ Xw (SR/600 >! ( SRI ) | (WF) | (FT) J (1> 

i * 

where SR Is the dally solar radiation, WF Is the ratio of dally 
evapotransplratlon to potential evapotransplratlon (Stuff et a|., 1978), 
and FT Is a dally temperature function that relates growth rate to 
temperature (Coelho and Dale, 1980). The term "solar radiation 
Intercepted" (SRI) refers to the component that relates the available 
energy to that Interacting with the leaves (Unvlll et a)., 1978). 

For the maize (corn) canopy, Unvlll et a|., (1978) has related SRI 
to leaf area Index (LAI) using fho following equation; 

SRI ° I - exp(-0.79 • LAI) (2) 

which has a range from zero to one and reaches 90 percent when LAI Is 
about 2.9, as shown In Figure 1. 

Investigators working In the optical region of the electromagnetic 
spectrum have related the ratio of an IR-band radiance and a red-band 
radiance to LAI to monitor plant performance for corn and wheat 
(Daughtry et al., 1980; Tucker et al., 1981). They have found that by 
Integrating this ratio over a given period of time, accumulated dry 
matter, which Is related to yield, may be estimated. 

A correlation of 0.89 was reported between LAI and "greenness" for 
corn by Daughtry et al., (1982), and a correlation of 0.84 between the 
true yield of corn and Its estimated yield, based upon meteorological 
data and spectral data, was also reported. 
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Figure 1. Solar radiation weighting factor for determining 
Interception of solar energy by crop as a 
function of its leaf-area index (from Daughtry 
et al . , 1982) , 
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3.0 VEGETATION CANOPY MODEL 


ORIGINAL PA® H 


In the microwave region, a vegetation canopy may be modeled as o 
lossy volume of scattering elements, bounded by air above and by a 
scattering soli surface below (Figure 2a). The backscatterlng 
coefficient o^an represents the sum of the contributions from the 

canopy Itself, from afreet backscatterlng by the soli (Including two-way 

* 

attenuation to account for propagation between the air-canopy boundary 
and the canopy-soil boundary and back), and from multiple scattering 
between the canopy scattering elements and the soli surface. 

In general, the scattering behavior of the canopy volume Is 
governed by the dielectric properties and geometric configurations of 
the scattering elements (leaves, stalks, and fruit), the latter 
being defined with respect to the wavelength, direction, and 
polari^f^loTi of the Incident wave. A first-order canopy backscatterlng 
mocH was developed by Attema and Ulaby (1978), who treated the canopy 
as a water cloud consisting of a collection of Identical water particles 
(Fig. 2b) character ized by a uniform scattering phase function. The 
water-cloud assumption Is a consequence of the domination of the 
dielectric constant of green vegetation (which Is a mixture of 
vegetative matter and water) by the dielectric constant of water; while 
the relative dielectric constant of water Is about 80 (below 10 GHz), 
the dielectric constant of dry vegetation is on the order of 2 to 3 
(Carlson, 1967; Jed I Icka et a|., 1983). Ignoring contr I but Ions 
resulting from multiple scattering between the canopy particles and the 
soil surface, the backscatterlng coefficient of the canopy (including 
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(T) Direct Backscatter from Plants 

(2) Direct Backscatter from Soil 

(Includes two-way attenuation by canopy) 



(a) Scattering Sources 


(D Backscatter from Vegetation "Water Particles" 



(b) Equivalent Cloud Model 


Figure 2. Backscattering contributions from a vegetation canopy, 
(a) scattering sources, and (b) equivalent "cloud" 
representation in terms of water scatterers. 
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soil contributions) Is given by: 

o 0 omcfNAL* pass rj 

can( 0 > * a veg<*) + ergo! I CO) OF POOR QUALITY ( 3 ) 

o o 

where a ve g is the contribution of the vegetation volume, a S oll I s the 

backscatter I ng contribution of the soil surface In the presence of 

vegetation cover, and 0 Is the angle of Incidence relative to naclr. 

Assuming the scattering water particles to be uniformly distributed 

within the canopy volume, Attema and Ulaby (1978) derived an expression 
o 

for a veg( Q ) by Integrating the backscatter I ng contributions of thin 
strata between the alr-vegetat|on boundary and the vegetation-sol 1 
boundary, 

o «r v cos 9 

%eg(® ) = Cl “ T 2 (Q)3# m 2 m” 2 

2 K e 

(4) 

where 

T 2 ( 6 ) = exp(-2x e h secO) (5) 

h = canopy height, m, 

<e = canopy extinction coefficient, m”l, 
a v = canopy volume backscatter I ng coefficient, m”l, 
and T 2 is the two-way transmission coefficient 
of the canopy, m 2 m“ 2 . 

The soli backscatter I ng coefficient, CT soi|(9)# depends on tne 
soil's surface roughness and its dielectric properties, the latter being 
governed strongly by the moisture content of the soil surface layer. 
When radar Is used as a vegetation monitor, the angle of incidence 0 
usually is chosen to be greater than 40° and the wavelength Is chosen to 
be on the order of 3 cm or shorter. Both choices are made, In part, in 
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order to make T as small as possible, thereby decreasing the soli 

backscatter Ing contribution (the second term In (3)) to a negligible 

level In comparison tc the vegetation backscatter Ing contribution. If 

the second term In (3) Is Indeed much smaller than the first. It may be 

o 

neglected. In general, however, a S oII rcay be described by the simple 
expression (Attema and Ulaby, 1978): 

a so||(e) = CC( 0 ) m s ] (6) 

where C(0) Is a constant for a given wavelength and polarization 
configuration, and m s Is the soil moisture content. In general, C(e) Is 
a function of the soli surface roughness, but If A < 3 cm, C(e) is 
approximately roughness- I ndependent over the range of random roughners 
usually encountered In the case of agricultural crops. 


4.0 .GlfllUd,J:1g.Cle.l in. Terms of. Canopy Water Content 

Attema and Ulaby (1978) assumed that the "equivalent" scattering 
water particles of the vegetation volume are all spherical In shape. 
Identical in size, and small relative to the wavelength A, In which 
case, 

g v = N (7) 

K e = N Q e (8) 

where N(m’ 5 ) Is the number of scattering particles per unit volume and 
^(m 2 ) and Qe^nr 2 ) are the backscatter Ing cross-section and extinction 
cross-section for y single particle, respectively. For an atmospheric 
water cloud, k q [ s directly proportional to the cloud's volumetric water 
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content m v (kg rrf 3 ). Hence, for the canopy "cloud," < Q may be expressed 
as 

*e = A | m v (9) 

where A| Is considered a constant at a given wavelength In general, 
A] may also be a function of the physical temperature of the canopy 
(through the temperature dependence of the dielectric constant of 
water), the wave pol ar Izatlon, and the shapes and sizes of the real 
scattering elements (leaves, stalks, and fruit). Approaches that 
Incorporate seme of the geometrical characteristics of the canopy are 
discussed In the next section. 

Similarly, 

a v = ^2 % (1 0 ) 

where A 2 Is a constant, and 

°v A a 

2 k 0 2 A 2 (11) 

Hence, the above ratio, which appears as the front part of the 

expression given by (4), is a constant. This Is a consequence of the 

assumption that all the equivalent water particles are spherical In 

shape and Identical In size. If the particles are spherical, but their 

sizes are distributed over a range of values, the relation given by (9) 

will continue to hold, but that given by (10) may not. For an 

atmospheric cloud, for example, CT V is proportional to m2 rather than 

v 

to m v . Hence, in general, the ratio on rhe right-hand side of (11) may 
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be written In the form: 
o 

V * 

a A4 mv 

2 K e (12) 

where x is a constant exponent. 

Upon Inserting (4), ( 6 ), (9), and (12) Into (3) and (5), the 

following expressions are obtained: 

°can(9) = a 4 my cose Cl ~t 2 ( 0 )] 

+ T 2 (0) C ( 0 ) m s (13) 

and 

t 2 ( 0 ) = exp(-2 A) m v h sece) (14) 

where A], A 4 , C, and x are constants for a given crop type and sensor 
configuration (wavelength, polarization, and incidence angle), m v and h 
are physical parameters of the vegetation canopy, and m s , Is the soil 
moisture content. As was mentioned earl ler, Attema and Ulaby (1978) 
assumed that the particles are all Identical in size, which corresponds 
to the case where x = 0. By regressing experimental data against a 
model of the form given In (13), they determined the values of A], A 4 , 

and C for several crops at each of several frequencies. In a more 
recent study, Hoekman et al. (1982) also found the cloud model to 
provide a good description of the temporal behavior of the 
backscatter ing coefficient for several crop types. For cereal grains, 
however, the model failed to predict the large changes in o° an that 
were observed In conjunction with the appearance of the heads. This led 
Hoekman et al . (1982) to subdivide the vegetation layer into two layers: 
a lower layer representing the stalks and leaves, and an upper layer 
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representing the ears or heads. Using this two- layer form of the mode!, 
they obtained good agreement between measured and model-predicted values 
of over the full growth cycle for each of eight crops. 


5.0 0£S.QR.1£T1QILQF. EXEERimiT 

In an effort to Investigate further the connection between the 
physical parameters of the canopy (m y , h, LAI, etc.) and the radar 
backscattering coefficient a can> field experiments were conducted at a 
test site near Manhattan, Kansas, during the growing seasons of 1979 and 
I960. The experiments were conducted by the Remote Sensing Laboratory 
at the University of Kansas in cooperation with the Evapotranspi ration 
Laboratory at Kansas State University, Manhattan. 

o 

A truck-mounted scatterometer was used to measure a C an a "f 
13.0, 17.0, and 35.6 p GHz (or equivalently, X = 3.5 cm, 2.3 cm, 1.76 cm, 
and 8.4 mm). A total of 14 fields were observed In 1979, 
six each of corn and sorghum, and two of wheat (see Table 1). Each observation 
sequence consisted of measurements at angles of Incidence 0 of 


30°, 50 o, an< j 

700 

for 

each 

of three linear polarization 

conf Igurations 

(HH, 

HV, 

VV). 

In 1980, the observations 


were limited to 0 = 30°, and the number of fields was reduced to six, 
thereby Increasing the number of temporal observations per field from an 
average of seven for 1979 to 23 for I960. 

In support of and contemporaneously with the radar observations, 
several plant- and sol I -properti es were measured, including plant 
height, the fresh weight W w and dry weight W ^ of individual plant parts 
(leaves, stalks, and fruit), plus the density of the whole plant, the 

11 
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LAI, soil moisture content m s , planting density N p (plants m“ 2 ), and 
stage of growth. The volumetric water content of the canopy can be 
determined from*. 

m v * <W w - W<j) Np/h (15) 

where W w and W<j are the average wet and dry weights of an entire plant, 

0 

Examples of the observed temporal patterns of o can are shown In 
Figures 3 to 5 for winter wheat, corn, and sorghum, all at 13.0 GHz, 9 ■ 
50°, and VV polarization conf Iguratlon. Also shown In these figures 
are plots of the green-leaf area Index, which show a fair degree of 
similarity to the plots of ^can# except for the period preceding 
harvest. It was also observed from the field data that the 

leaf-area- Index (LAI) was highly correlated to the dry and to the wet 
(leaf) vegetation biomasses. This proves to be an Important link, as 
will become apparent In the next section. 


6.0 Multi-Component Cloud Model 

After completion of the processing and calibration of the data 
described In the previous section, the data were used to evaluate the 

applicability of the model given by (13). First, the model was 

0 

evaluated by correlating the measured value of a caq to the value 

calculated using (13) with the constants Aj, A 4 , and C being assigned 
the values that were determined by Attema and UJaby (1978). This 
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WHEAT FIELD #1 1979 13.0 GHz 50 Deg VV 



14 


JULIAN DATE (1979) 

figure 3. Measured temporal patterns of the backscattering coefficient and leaf-area index (LAl) 
of a wheat field. Stage of growth is indicated by the Feekes scale (Large, 1954} . 
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JULIAM DATE (1980) 

Figure 5. Measured temporal patterns of the backscattering coefficient and leaf-area index (LAl) 

of a sorghum field. Stage of growth is indicated by the Vanderlip growth stage (Vanderlip, 1972) 
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resulted Iri a linear correlation coefficient p between 0.6 and 0.8 (for 

the various combinations of crop type, microwave frequency, angle of 

Incidence, and polarization). Next, another correlation analysis 

between measured and calculated values of c£an was performed, except 

this time the constants were first determined empirically by regressing 

the data against the model. Some Improvement In the magnitude of p was 

obtained, but p exceeded 0,8 only In a few cases. 

To Improve the applicability of the cloud model, two modifications 

were made, both of which resulted In better agreement between the 

0 

measured temporal pattern of a can and that calculated by the model. 


6.1 Corn and Sorghum 

The first modification provides for separate accounting for the 
radar backscatterlng contributions of the leaves, the stalks, 
and the soil. For corn and sorghum, the canopy was assumed to consist 
of two layers: an upper layer of height hj, dominated by leaves, and a 

lower layer of height h£, dominated by stalks (Figure 6a). The 
backscatter contribution of the fruit was Ignored, In part to simp I Ify 
t (. i model, and In part because of the results of a defoliation 
experiment for corn which showed that the backscatterlng contribution of 
the fruits Is much smaller than that of the stalks or leaves (Ulaby, 
1982). These assumptions lead to: 

^can^S^ = a |eaf (16) 

+ CT stalk< 9 > * °so! I 
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where ^ | ©af (9) has the same form as (4), OF POOR QUALITY 

A 0 . COS0 5 

n - T|Mf (0)l, <17. 

2 

and T |eaf(S) * exp(-2 k 0ji hj sect)). (18) 

in the above expressions, hj Is the "effective" height of the top 
layer, assumed to consist of leaves exclusively, and 0 VJ< and K e & are the 
volume backscatterlng coefficient and extinction coefficient of that 
layer. The stalk's contribution was assumed to be given by: 

a stalk^^ n Tfeaf^ 9 ) °stalk( 0 l hi 8 0) 

(19) 

where ^sfaik (hi B 0) Is the stalk contribution In the absence of 
leaves. Finally, the soli component is given by a form similar to that 
of (6) except that the ‘attenuation through the canopy is now due to two 
layers: a top layer with transmission coefficient T | e af and a lower 

layer with transmission coefficient T s + a ||<, 

4 

CT sol|(9) * Cc(9) m s 2 T |e a f ( 0) T 2 ( 9 ) . 

stalk' ' f (20) 

where 

2 

T stalk< 6) 3 exp <“2 k 6s h 2 sec9) (21) 

and k qs |s extinction coefficient of the stalk layer. 

The second modification to the model involves the use of green LAI. 
0 

In Section 4, a can was expressed In terms of three physical parameters: 
the soil moisture content m s , the volumetric wafer content of the canopy 
m v , ana the canopy height h. Actually, the key quantity driving the 
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model Is the product m v h, which represents the total amount of 

Integrated water contained In a vertical column of unity horizontal 
cross-section. The model In this section divides the canopy Into two 

layers, with the top layer assumed to contain essentially all the leaves 
and the bottcm layer turned to contain the entire length of the stalk. 
The quantity myhi for the top layer Is related to the wet and dry 
biomasses of the leaves, which In turn, are related to the green LAI 

(Zrust et al., 1974; Ashley et at., 1965; Aase, 1978; Holben at al., 

1980). Hence, Instead of expressing T| e af In terms of myhi (as In 
(14)), It will be defined In terms of LAI: 

Tjeafte) * exp (-2 B tea f LAI sec6). (22) 

For the stalk layer, it will be assumed thaT the scattering by the 
stalks, as well as the attenuation by them, is proportional to m v h 2 
where m v Is the volumetric water content of the stalk layer, hence, 
o 2 

a stalk(®) = A s fal(< m v h 2 T| ea f(6) sine (23) 

and 

2 

T stalk< 0 ) = exp(-2 Q s -ta i k %h2 secs) (24) 

where Asfgj^ and B s t a ||< are constants (for a given crop type, microwave 
frequency, and polarization). Finally, It will be assumed that the 
ratio a v Q heys the assumptions that led earlier to (11). Grouping 
the above terms leads to*. 

a can< 9 > 58 A| eaf cos0Cl - t2 / Q \ i 

leaf ' ' J 

Astalk m v h 2 ? f ea f( 0 ) sin 9 

* Ksoil< 9 > "Lf (9 > T s t alk (l3) <25) 
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The above expression contains five coefficients and three physical 

parameters: (a) green LA I, (b) m v h 2 , the total water content of the 

stalks per unit area, and (c) m s , the volumetric soil moisture content. 

The coefficients were determined by using a non-l Inear regression 

program that minimizes the least-squares error between the predicted and 

observed values. Data from all three fields (of 1980 sorghum or 1980 

corn) were used to generate the coefficients, and then the model was 

used to predict the temporal pattern of a 0 for each field Individual |y« 

The values of the coeffclents obtained are given In Tables 2 and 3 for 

corn and sorghum, respectively. Also given Is the linear correlation 

o 

coefficient between the measured and predicted values of ° can for each 

field at each of the four frequencies. Figures 7 and 8 show comparisons 

for corn and sorghum, respectively. Each figure contains a plot of the 
o o 

predicted a can/ denoted 0 p r ed> as wel I as plots of Its component 

contr lout Ions as defined by (16) or by the three individual terms In 
(25). It Is observed that: 

(a) Except for the period prior to Julian day 170 (for which the 
pi ant height was less than 0.8 m for corn and less than 0.6 m for sorghum) 
and for the pre-harvest period after Julian day 240, the leaves provide 
the overwhelming majority of the backscattered energy. 

(b) The soil and stalk terms are Important only if LAI < 0.5. 

(c) The stalk attenuation coefficient, B s f a ||<, was found to be very 
small In magnitude, and when It was set equal to zero (I.e, setting 
^stalk = 1)» the model remained effectively the same. 

(d) A better fit could be obtained for corn if the leaf attenuation 
coefficients, B| ea f> were to be assigned different values for the two 
periods before and after Julian day 190, which corresponds to the date 
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TABLE 2 

Coefficients for the Corn Model Given by Equation (25) Rewritten in 
the Following Abbreviated Form for 0 = 50° and VV Polarization: 




vC 

e'- 

en 

a II 




li- 

03 






re 

•M 




— — 


b 



, — . 


t— 


C4 — 

iH 


CM 


0 





x 

4-1 

3 « 



• 



x 

b 


■ 

4- 


• 

01 


r™ * 


03 

. — . 


•— 

+ 

4- 

CM 

O 


> 

CD 


03 

x 

Oi p— 

< 

E 

x 

JZ 

<U 



_l 



r— 

CN *— 

• 



• 

>N 

re 

e- 


• 

• 


a. 

44 


> 



.X 

0 

3 10 

1 

e 


ti- 

— 

c 

O 



E 

re 

re 

re 


r** 

* 


ai 

4J 

O 

+ 

i — i 


• 


- to 




_x 


m 

03 

b 

14- 

4* 



1 

1 

x 

re 

03 

re 


'• — ' 


4J 

a) 

<D 

4-1 

0 

a 

a 




- to 

to 

X 

X 


b 

< 

< 

O 

b 

b 

X 


II 

c 

m 

o o 

o 


ii 

ti- 

ro 

<u 

CVJ — 
6 - 

b 

U 


JX 

re 

4-> 

W l/| 
£-* 


— CN 


*o 

c 

re 


a ) 

x 

3 



(0 

• 

X 

E 

0 0 

re 

b 

u 

• 

Ol 

to "o 

0 

b b 

L- 

3 C 

a 

r— •- 


re x 

c 

> E 

0 

0 

»— » 

TJ u 

to 

b 

10 

> to 

b 

C -0 

u 

b r- 

b) 

to b 

b 

x ■— 

u 

O 4— 

u 

b b 

re 

X b 

b 

4J L. 

c 

X 

• — 

C 4-1 

— 

b 

c 

4) »■— • 

0 


c 

4-» re 


b 

re 

X t~ 


0 

4 - 

C ti- 

0 

0 

— to 

b 

4J — 

to 

re 

3 

0mm 


b — 

b 

L. 4— 

X 

u re 

4-1 

o a 
u 

r~ 

■b 

CJI 

L- c 

3 

re re 

O 

b 

1_ 

c - 

X 

«— ■ — 

4-4 

— — 


•b 

"b 

b — 

b 

X b 

C 

4J •— 


U- 

*i 

(0 

L. 

— u 

b 

0 

4J 

— 4- 

b 

a. *b 

■b 

b 


4-1 U 

b 

C O CL 

1— 

b b 

b 

I— 

■— 

u to 
— a) 

to 

ti— 3 

4-1 

ti_ — 

c 

b ra 

b 

0 > 

• — 

u 

0 

■b 

*— 

C b 

14 - 

O 4-1 

w- 

O 

b 

4-4 •— 

O 

10 TJ 

y 

— 0) 


b t- 

b 

1 - Q. 

X 

i_ 


O b 
U X 


4-1 


b 

X ~c 

b- c 


re 





CM 


cN 

CM 


a 

II 


2 




> 


4— 


CM 



a 

II 


2 






tn 

— 

M3 

r— 


re 

II 

a 



•w* 


H 

r— 

1 

03 

*5 

m 

to 

E 

0 

u 


■H 

■X 

1 

— 

cn 

re 


4-4 


- w 

<N 

03 

E 

ti- 


re 


b 


— — 


cn 


-X 


— 


re 


4-1 


- 10 


<c 


ti- 


ro 


b 




< 


>• 


u 


c 


b 

N 

3 

X 

era 1 

b 


i_ 


u. 



22 


8.6 0.1359 0.01662 1.046 0 0.2118 0.895 0.837 0.931 0.895 

13.0 0.1697 0.01783 1.1 2^ 0 0.2094 0.885 0.900 0.899 0.928 

17.0 0.1925 0.0125** 0.895 0 0.271 0.852 0.845 0.860 0.938 

35-6 0.2209 0.02487 0 . 8 *i 30 0 0.1**51 0.914 0.894 0.938 0.926 


TABLE 3 

Coefficients for the Sorghum Model Given by Equation (25) Rewritten t 
the Following Abbreviated Form for 0 = 50° and VV Polarization: 
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at which the LAI reaches Its peak value. It has been suggested that the 
optical transmission properties of leaves change In magnitude as the 
plant enters the senescence phase. The same may be true for the 
microwave part of the spectrum, but no data exists at present to support 
such a hypothesis. 

6.2 MhflaJ: 

Wheat has a significantly different geometry from that of corn or 
sorghum, thereby necessitating that certain modifications be Introduced 
to the model of the previous section. First, the wheat stalks are much 
smaller and contain a smaller portion of the total plant water. Hence, 
the stalk's backscatter Ing contr Ibutlon w 1 1 I be neglected. Secondly, 
the size, location, and relative water content of the heads suggest that 
their backscatter Ing contributions should be accounted for explicitly 
(rather than being Ignored or lumped with the stalks, as was done In 
the previous section for corn and sorghum). Additlonaly, since the 
heads of the wheat plants are above the leaves (Figure 6b), the leaves 
should not attenuate the backscatter Ing from the head; rather, It should 
be the other way around. These modifications, plus seme others 
discussed below, lead to the expression 

C can* 9 * ~ a leaf< a > + CT head< 0 ) + so 1 1 ( 0 ) 

- *|«f ' LA| C ‘ ’leaf 1 "’ 3 T head 9 
+ A head (0) ' M head 
♦ c 5 o ,,(0> m s T 2 eaf ( e ) T 2 eaa (0) 
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whore T| ea f( f )) Is given by (22), M head is the head biomass In kg m“ 2 , 

and 

T head <8> x ex P<- 2 B head M head sec8) - <27) 

The coefficient of the oJ ea ^ term Is (A )ea f • LAI), which Is somewhat 
equivalent to (12) with x » 1 . This form Is used, rather than simply 
A| e af» because It resulted In significantly better agreement between the 
measured and predicted values of o° an for wheat (as discussed In more 
detail In the next section). 

The field measurements did not Include direct measurements of the 
head’s biomass, M^ead? therefore, an estimate was needed. If the plant 
is assumed to be fully developed before heading takes place, then the 
dry weight due to the leaves and stalks should remain constant for the 
period during and after heading (assuming no leaves fall off the plant). 
Any change In dry weight after the onset of heading will therefore be 
directly related to the development of the head. This reasoning leads 
to: 

* 

M head(t) = \ 0, for t < to; to = heading date 

|wq(t) - W<j(t 0 ), for t > t 0 (28) 

where is the total plant dry biomass In Kg m“ 2 . For the fields 
observed in 1979, t Q = 136 (May 16, 1979). 

Following a determination of the values of the five coefficients 
employed In the model (see Table 4), plots were generated to compare 
°pred +0 a obs* E * am P ,es are shown in Figures 9a and 9b for 13 GHz and 
35.6 GHz, respectively. In both cases, the model provides a good fit to 
the measured data, and Indicates that the Dackscafter Is almost 
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TABLE k 

Coefficients for the Wheat Mode) Given by Eq. (26), Rewritten in the Following 
Abbreviated Form for 8 = 50° and VV Polarisation: 
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The coefficients were determined through the use of a nonlinear regression program. 
The correlation coefficient pj defines the linear correlation between the observed 
values a 0 , and the predicted values o° , for Field i. 

nhc * 


■ F1RD #1 13.0 GHz W 50 Deg 
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OF POOR QUALITY 



Figure 9. Comparison of measured backscattering coefficient of wheat with that computed 
using Eq. ( 26 ) at (a) 13-0 GHz and (b) 35.6 GHz, 50° incidence angle, and VV 
polarization. 




WHEAT FIELD #1 35.6 GHz W 50 Deg 
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exclusively due to the leaves between Julian days 120 and 155. After 
Julian day 155, green LAI goes to zero, and the observed backscatter 
becomes due entirely to the heads and the soil surface. 

7.0 APfftQXl M ATEJWl. t N,T£R M £-QF— LA 1 

The models discussed In the previous section appear to agree well 
with experimental observations, which makes them useful tools for 
evaluating the sensjtlvly of the radar backscatter I ng coefficient to the 
physical parameters of the plant canopy and soil surface. From an 
applications standpoint, however, radar would be far more useful If the 
backscatter I ng coefficient It measures could be used to estimate a 
single, but Important, physical parameter of the canopy, such as green 
LAI. The conclusions arrived at In the previous section relevant to the 
temporal "signature 1 ' of °° can way be summarized as follows: 

(a) During the early stage of growth (short plants and LAI < 0.5), 

the magnitude of a 0 may be affected, even dominated, by sol l-molsture 
can 

conditions. 

(b) During the stage of growth character Ized by LAI > 0.5, Is 

can 

dominated by the leaf contributions. Hence, cr° s cr? - for 

' can leaf 

LAI > 0.5. 

(c) During the stage of growth prior to harvest and LAI < 0.5, 
o 

a can Is dominated by the soil and stalk contrl but Ions for corn and 
sorghum, and by the soil and head contributions for wheat. 

Based on the above conclusions, the models may be simplified as 
presented In the next sections. 
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7,1 .CQ r . Q ,^ntLSacshum 

Ignoring the presence of the stalks In the canopy, (22) and (25) 
may be combined and simplified to: 


o < I 

a can( 9 * * A| eaf + S(6) exp(~ B|eaf 


LAI) 


( 29 ) 


where: 

8 leaf 2 B | ea f secB (30) 

i 

S(6)» j C sot | (0) m s - A | eaf » for LAI < 0-5 (31) 

| C soil (e> ■ A iaaf . f °i-LAI>0.5 (32) 

In (32), Csoii = C so j|*m s f where m g I s some representative val ue 

of the mean soli moisture content. It may be taken as 

m s * 0.20 g cm 3. For LAI > 0.5, the soli contribution may be assumed 
to be a constant (rather than Ignored), thereby simplifying o° an + 0 
the point where It becomes a function of only one variable: green LAI. 

t Q 

The effect of soil moisture variations on n can [ s illustrated In Figure 

o 

10, which shows the range of values that a can may take In response to a 
change In soil moisture content from 0.05 g cm' 3 , representing 
relatively dry conditions, to 0.25 g cm" 3 , representing relatively wet 
conditions. Clearly, soil moisture variations ar, important only If LAI 
I s smal I . 

Figures It and 12 present experimental observations of for 

corn and sorghum, respectively, plotted as a function of LAI. In each 
case, the data Includes all observations made over the season for all 

t 

three fields of that crop type. Also shown are plots of the 
single-parameter (LAI) model (as defined by (29) and (32)) fitted to the 
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data. The observed scattering Is due to several sources, Including; 
(a) sensor measurement precision, estimated to bo about ±1 dB, (b) soil 
moisture variations, <c) errors associated with the measurement of LAI, 
and (d) wlthln-fleld and between-f lei d variations that are not accounted 
for by the model or the data. 

7.2 m tal 

For wheat, a si ngl e-parameter model (In terms of LAI) may be 
written In the form; 

a can<9> = A|eaf^AI) x Cl - exp(-B, l eaf . LAD] 

+ C so j|( 0 ) exp(- B| ea f • LAI), ( 33 ) 

for LAI > 0.5 

which reduces to the combination of (29) and (3D If x = 0. Actually, 
this form was used In the non-linear regression program for all three 
crops, but the exponent, x, was assigned a value of zero or very close 
to zero by the program In each and every case for al I the corn and 
sorghum data sets. For wheat, on the other hand, x wa a assigned a value 

between 0.7 and 1.4, depending on the frequency and polarization 
combination. The goodness of fit of (33) to the measured data Is 
i I I ustrated in Fig. 13. 

8.0 QQNCLUS.IQNS 

By extending the canopy cloud model of Attema and Ulaby (1978) from 
Its Initial consideration of the canopy as a single layer with uniform 
volume scattering and attenuation properties to one consisting of more 
than one layer, it was possible to examine the relative backscatTer i ng 
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contributions of the leaves, stalks, and heads (of wheat). The analysis 
of radar backscatter Jng data measured at an Incidence angle of 50° and 
microwave frequencies of 8.6, 13.0, 17.0, and 35.6 GHz Indicates that 
the backscatterlng coefficient of the canopy Is dominated by the leaf 
contribution If green LAI Is greater than approximately 0.5 for corn and 
sorghum and If, additionally, the heads have not yet appeared (for 
wheat). During the early stage of growth, a soli backscatterlng 
contribution may be very Important, and for the end-period prior to 
harvest, the backscatterlng contributions of the soli and the stalks are 
Important for sorghum and corn, and of the heads and soil for wheat. A 
simplified version of the model was developed In terms of a single 
parameter: green LAl. 


40 


ACKNOWLEDGMENT 


This research was supported by the National Aeronautics and Space 
Administration, Johnson Space Center, under Contract NAS 9-15421. 


REFERENCES 


Aase, J. K. (1978), Relationship between leaf area and dry matter 
In winter wheat, Agron... J.. 70:563-565. 

Ashley, 0. A., B. D. Doss, and 0. L. Bennett (1965), Relation of 
Cotton leaf area Index to plant growth and fruiting. Agron. J. 
57:61-64. 

Attema, E. P. W., and F. T. Ulaby (1978), Vegetation modeled as a water 
cloud. Radio Science . 13(2) :357-364. 

Brisco, B., and R. Protz (1980), Corn field Identification accuracy 
using airborne radar Imagery. Can. ... I. -Remote- Se.ns.ljip . 6(1) : 1 5-25 . 

Bush, T. F., and F. T. Ulaby (1978), An evaluation of radar as a crop 
classifier. Remote -Sensing, ..EnYlcfiri,. 7:15-36. 

Carlson, N. L. (1967), Dielectric constant of vegetation at 8.5 GHz. 
Ohio State University Electroscience Laboratory Technical Report 
1903-5. Columbus, Ohio. 

Coe I ho, D. T. , and R. F. Dale (1980), An energy-crop growth variable and 
temperature function for predicting corn growth and development: 
planting to silking. Agron. J. 72:503-510. 

Dale, R. F. (1977), An energy-crop growth variable for Identifying 
weather effects upon maize yields. In Agrometeoro I ogy. of Maize 
(Corn) -.Crop, W. Baler, R. H. Shaw, L. M. Thomson, and 
R. E. Fetch (eds.). fr.oc . ,-Sy.mp ...Aflr.Qine.t9QCfi L-Ma-tee. .... S Corn) 

Crop . World Meteoroi. Org., Iowa State Unlv., Ames, Iowa, 
5-9 July 1976, WM0 No. 481, Geneva, Switzerland. 

Daughtry, C. S. T. , and N. C. Fuhs (1980), Determination of the value of 
spectral Information In estimation of agronomic variables 
associated with yields of corn and soybeans. A gRl STARS Report 

Daughtry, C. S. T., K. P. Gallo, and M. E. Bauer (1982), Estimation of 
agronomic variables associated with yields of corn and soybeans. 

A g R 1 S T ARS - Rep a r.t. , .SRr P2rQ4266 , pp. 31-42. 

Hanway, J. J. (1971), How a corn plant develops. Special Report 
No. 48, Iowa State Unlv., Cooperative Extension Service, 

Ames, Iowa. 

Hoekman, D. H. , L. Krui, and E. P. W. Attema (1982), A multilayer model 
for radar backscatterlng from vegetation canopies. Digest. 2nd 
Annual ialL — Sao.scl.ence and Remote Sensing Svmp. r Vol . II, Munich, 
West Germany, 1-4 June (IEEE Cat. No. 82CH14723-6) . 


Hoi bon, B. N. , and C, Fan (1980), Spectra) assessment of soybean 

leaf area and leaf biomass. Photoqramm. Eng, and Remote Sensing . 
46(5) :651-656. 

Hoogeboom, P. (1982), Classification of agricultural crops In radar 

Images. Digest. - 2nd Annu al Inti. Geoscience and Bfimsla 

Sensing Symp. r Vo! . II, Munich, West Germany, 
1-4 June (IEEE Cat. No. 82CHI 4723-6) . 

Jedl Icka, R. P., and F. T. Ulaby (1983), Measurements of the radar 
microwave dielectric and attenuation properties of plants. 

Pros .-..Nation a I Radio Sc lanes,.,. Mooting ILlRSi GamuJssIan — El, 

5-7 Jan., Unlv. Colorado, Boulder, Colorado. 

Large, E. C. (1954), Growth stages In cereals, 1 1 1 ustratlon of the 
Feekes scale. Plant Pathol.. 3:128-129. 

LI, R. Y., F. T. Ulaby, and J. R. Eyton (1980), Crop classification with 
a Landsat/Radar sensor combination. Presented at the 6th Purdue 
Symp. on Machine Processing of Remotely Sensed Data. 2-6 June, 

Purdue Unlv., West Lafayette, Indiana. Supported by 
NASA Contract NAS 9-15421. 

Llnvllle, D. E., R. F. Dale, and H. F. Hodges (1978), Solar radiation 
weighting for weather and corn growth models. Agron. J. 

70:257-263. 

MacOonald, R. B., and F. G. Hall (1980), Global crop forecasting. 
Science . 208:670-679. 

Shanmugan, K. S., F. T. Ulaby, V. Narayanan, and M. C. Dobson (1983), 
Identl f fcatlon of corn fields using multidate radar data. Remote 
Sens.Inq of. Environment. Forthcom I ng. 

Stuff, R. G., and R. F. Dale (1978), A soli moisture budget model for 
accounting for shallow water table Influences. 
SqI l....Sci^.Soc.,-Am.„.J. 42:637-643. 

Tucker, C. J., B. N. Hoi ben, J, H. Elgin, Jr., and J. E. McMurtrey III 
(1981), Remote sensing of total dry matter accumulating In winter 
wheat, Remote. -Sens-Lnc of Environment. 11:171-189. 

Ulaby, F. T. (1982), Review of approaches to the Investigation of the 
scattering properties of material media. Digest. 2nd Annual Inti. 

Geoscience and Remote S ensing Svmp. Munich, West Germany, 

1-4 June, (IEEE Cat. No. 82CH1 4723-6) . 


Vanderj Ip, R. L. (1972), How a sorghum plant develops. Contribution 
No. 1203, Agron. Dept., Kansas Agricultural Experiment Station, 
Manhattan, Kansas. 

Zadoks, J. C. , T. T. Chang, and C. F. Konzak (1974), A decimal code 
for the growth stages of cereals. Weed Research . 14:415-421. 

Zrust, J.,E. Partykova, and J. Necas (1974), Relationship of leaf area 
to leaf wleght and length In potato plants. Photosynthet lea 
(Prague). 8:118-124. 


APPENDIX A 
DATA DOCUMENTATION 

Observed values of the backscatterlng coefficients and ground-truth values. 


VARIABLE 

UNITS 

DESCRIPTION 

SIGMA 0 

dB 

backscatterlng coefficient 
measured, pol * VV, 9 «* 50 ° 
(missing value represented by 0) 

Plant Height ■ 

m 

canopy height 

Plant Water 

-3 

kg m 

volumetric normalized plant 
water content 

Soil Moisture 

g cm 

volumetric soil moisture 

Leaf Area Index 

2 -2 
m m 

green leaf area index 

Plant Dry Mass 

-2 

kg m 

total canopy dry biomass per 
square meter 

Leaf Wet Mass 

-2 

kg m 

fresh mass of leaves per square 
meter 

Leaf Dry Mass 

-2 

kg m 

dried mass of leaves per square 
meter 

Stalk Wet Mass 

kg m" 2 

fresh mass of stalks per square 
meter 

Stalk Dry Mass 

kg m 2 

dried mass of stalks per square 
meter 
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